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Description 

[INTEGRATED CIRCUIT YIELD 
ENHANCEMENT USING VORONOI 

DIAGRAMS] 

Cross Reference to Related Applications 

[0001] The present application is related to pending U.S. Patent 

Application 10/ , filed concurrently herewith to Allen 

et al., entitled "CRITICAL AREA COMPOSITE FAULT MECH- 
ANISMS USING VORONOI DIAGRAMS" (IBM Docket No. 
BUR920030136US1). The foregoing application is as- 
signed to the present assignee, and is incorporated herein 

by reference. 
Background of Invention 

[0002] Field of the Invention 

[0003] The present invention generally relates to measuring criti- 
cal area in integrated circuit design, and more particularly 
to a method that uses Voronoi diagrams to measure criti- 
cal area as the design layout is changed. 



[0004] Description of the Related Art 

[0005] within this application several publications are referenced 
by Arabic numerals within parentheses. Full citations for 
these, and other, publications may be found at the end of 
the specification immediately preceding the claims. The 
disclosures of all these publications in their entireties are 
hereby expressly incorporated by reference into the 
present application for the purposes of indicating the 
background of the present invention and illustrating the 
state of the art. 

[0006] Advanced deep sub-micron technology enables millions 
of transistors to be fabricated on a single die. While this 
capability grants performance, the smaller size and higher 
density of layout features adversely affect yield [See: Pa- 
padopoulou, E., "Critical area computation for missing 
material defects in VLSI circuits ," Computer-Aided Design of 
Integrated Circuits and Systems, IEEE Transactions on , Vol. 20, 
No. 5, pp 583-597, May 2001; Fook-Luen Heng and Zhan 
Chen. "VLSI Yield Enhancement Techniques Through Lay- 
out Modification. "IBM T. J. Watson Research Center; and A. 
Venkataraman and I. Koren. "Trade-offs between Yield 
and Reliability Enhancement."Proc. of the 1996 IEEE National 
Symposium on Defect and Fault Tolerance in VLSI Systems, pp. 



67-75, November 1996]. Both the manufacturing process 
and geometry of the layout contribute to this loss of yield. 

[0007] a part of this yield loss comes from random defects. 
These defects occur during the manufacturing process 
and cause electrical faults when the chip is active. The 
types of electrical faults that may result include, but are 
not limited to, short-circuits, wire-breaks, and via ob- 
structions. The probability of yield loss due to random 
defects relates to the critical area of the layout, which can 
be measured using a Voronoi diagram technique [See: Pa- 
padopoulou, E. and Lee, D.T., "Critical area computation 
via Voronoi diagrams "Computer-Aided Design of Integrated 
Circuits and Systems, IEEE Transactions on , Vol. 18, No. 4, pp 
463-474, April 1999]. 

[0008] Critical area of a very large scale integration (VLSI) layout 
is a measure that reflects the sensitivity of the layout to 
defects occurring during the manufacturing process. Crit- 
ical area is widely used to predict the yield of a VLSI chip. 
Yield prediction is essential in today's VLSI manufacturing 
due to the growing need to control cost. Models for yield 
estimation are based on the concept of critical area which 
represents the main computational problem in the analy- 
sis of yield loss due to random (spot) defects during fabri- 



cation. Spot defects are caused by particles such as dust 
and other contaminants in materials and equipment and 
are classified into two types: "extra material" defects 
causing shorts between different conducting regions and 
"missing material" defects causing open circuits. 

[0009] | n some defect modeling techniques, defects are modeled, 
consistently, as circles. The underlying reason for model- 
ing defects as circles is the common use of Euclidean ge- 
ometry. The distance between two points, usually, is mea- 
sured by the length of the line segment joining the two 
points. This is the Euclidean distance. The locus of points 
a unit distance from a center point is usually called the 
"unit circle". In Euclidean geometry, the "unit circle" is a 
circle of radius one. 

[0010] | n reality, spot defects are not necessarily circular. They 
can have any kind of shape. Therefore, it seems appropri- 
ate to use other geometries if the critical area computa- 
tion can be simplified by modeling defects as squares, di- 
amonds or octagons. For practical purposes, a circular de- 
fect can certainly be approximated by a regular octagon. 
Yield estimation should not considerably depend on which 
of the above geometries is used to model defects as long 
as the geometry is chosen consistently. Therefore, the ge- 



ometry used for a particular computation, preferably, 
should allow critical area computation in the most effi- 
cient way. 

[001 1] a Voronoi diagram can also be used to enhance the com- 
putation of critical area. A Voronoi diagram of a set of 2D 
geometric elements (polygons, line segments, points) is a 
partition of the plane into regions representing those 
points in the plane closest to a particular geometric ele- 
ment. Here, "closest"is defined in terms of an appropriate 
geometry as mentioned above. These regions are called 
Voronoi cells, each of which is associated with its defining 
geometric element, called the owner of the cell. The set of 
points which separates two Voronoi cells is called a 
Voronoi bisector. The point where three or more Voronoi 
bisectors (or Voronoi cells) meet is called a Voronoi ver- 
tex. 

[0012] Based on the circuit design and under an appropriate ge- 
ometry, Voronoi diagrams can be constructed to model 
the effect of extra-material and missing-material spot de- 
fects. The Voronoi diagram partitions the circuit design 
into Voronoi cells within which defects that occur cause 
electrical faults between the same two shape edges in the 
design. This information can then be used to compute 



critical area, (e.g., see U.S. Patents 6,317,859, 6,247,853, 
and 6,178,539, which are incorporated herein by refer- 
ence). 
Summary of Invention 

[0013] The invention provides a method of calculating critical 

area in an integrated circuit design. Starting with an initial 
integrated circuit design, the invention associates vari- 
ables with the positions of the edges in the design. The 
invention associates cost functions involving the variables 
with the spacing among (between) these edges. The cost 
functions are in terms of critical area contributions. Criti- 
cal area contributions comprise a measure of electrical 
fault characteristics of the spacing among edges. The in- 
vention optimizes the position and length of the edges to 
reduce critical area contribution cost in a first direction 
across the integrated circuit design to produce a revised 
integrated circuit design. Then, the invention optionally 
repeats this process with the revised integrated circuit de- 
sign in a second direction to further reduce critical area 
contribution cost. 

[0014] The process of associating cost functions maps points 

within the spacing among edges in the design to the size 
of defects at those points that trigger an electrical fault, 



forming Voronoi cells. These Voronoi cells define Voronoi 
bisectors and Voronoi vertices which encode the variables 
associated with the edges. The invention defines cost 
functions based on these Voronoi elements but indepen- 
dent of their geometry. The cost function models critical 
area contributions of the edges in the design layout as 
they change position and length in a continuous manner. 
[0015] These, and other, aspects and objects of the present in- 
vention will be better appreciated and understood when 
considered in conjunction with the following description 
and the accompanying drawings. It should be understood, 
however, that the following description, while indicating 
preferred embodiments of the present invention and nu- 
merous specific details thereof, is given byway of illustra- 
tion and not of limitation. Many changes and modifica- 
tions may be made within the scope of the present inven- 
tion without departing from the spirit thereof, and the in- 
vention includes all such modifications. 
Brief Description of Drawings 

[0016] The invention will be better understood from the following 
detailed description with reference to the drawings, in 
which: 

[0017] Figures 1A and IB illustrate edges of design shapes and 



he Voronoi bisectors that are created between the design 
ihapes; 

: igure 2 is a Voronoi cell in three dimensions; 

: igure 3 is a diagram showing Voronoi bisectors; 

: igure 4 is a diagram showing Voronoi bisectors; 

: igure 5 is a diagram showing Voronoi bisectors; 

: igures 6A-6C are diagrams showing Voronoi bisectors; 

: igures 7A-7G are diagrams showing Voronoi bisectors; 

: igure 8 is a flow diagram illustrating processing accord- 
ng to the invention; 

: igures 9A-9C illustrate bisectors, and graphical cost 

unctions associated with bisectors; 

: igures 10A-10C illustrate bisectors, and graphical cost 

unctions associated with bisectors; 

: igures 11A-11C illustrate bisectors, and graphical cost 

unctions associated with bisectors; 

: igures 12A-12C illustrate bisectors, and graphical cost 

unctions associated with bisectors; 

: igures 13A-13C illustrate bisectors, and graphical cost 

unctions associated with bisectors; 

: igures 14A-14C illustrate bisectors, and graphical cost 



functions associated with bisectors; 
[0031] Figures 15A-15C illustrate bisectors, 

functions associated with bisectors; 
[0032] Figures 16A-16C illustrate bisectors, 

functions associated with bisectors; 
[0033] Figures 17A-17C illustrate bisectors, 

functions associated with bisectors; 
[0034] Figures 18A-18C illustrate bisectors, 

functions associated with bisectors; 
[0035] Figures 19A-19C illustrate bisectors, 

functions associated with bisectors; 
[0036] Figures 20A-20C illustrate bisectors, 

functions associated with bisectors; 
[0037] Figures 21A-21C illustrate bisectors, 

functions associated with bisectors; 
[0038] Figures 22A-22C illustrate bisectors, 

functions associated with bisectors; 
[0039] Figures 23A-23C illustrate bisectors, 

functions associated with bisectors; 
[0040] Figures 24A-24C illustrate bisectors, 

functions associated with bisectors; 
[0041] Figures 25A-25C illustrate bisectors, 

functions associated with bisectors; 



and graphical cost 



and graphical cost 



and graphical cost 



and graphical cost 



and graphical cost 



and graphical cost 



and graphical cost 



and graphical cost 



and graphical cost 



and graphical cost 



and graphical cost 



[0042] 


Figui 


-es 26A-26C illustrate bisectors, and graphical cost 




fund 


tions associated with bisectors; 


[0043] 


Figui 


es 27A-27C illustrate bisectors, and graphical cost 




fund 


tions associated with bisectors; 


[0044] 


Figui 


'es 28A-28C illustrate bisectors, and graphical cost 




fund 


tions associated with bisectors; 


[0045] 


Figui 


-es 29A-29C illustrate bisectors, and graphical cost 




fund 


tions associated with bisectors; 


[0046] 


Figui 


-es 30A-30C illustrate bisectors, and graphical cost 




fund 


tions associated with bisectors; 


[0047] 


Figui 


-es 31A-31C illustrate bisectors, and graphical cost 




fund 


tions associated with bisectors; 


[0048] 


Figui 


"es 32A-32C illustrate bisectors, and graphical cost 




fund 


tions associated with bisectors; 


[0049] 


Figui 


-e 33 illustrates Voronoi bisectors; 


[0050] 


Figui 


'e 34 illustrates Voronoi bisectors; 


[0051] 


Figui 


e 35 illustrates Voronoi bisectors; and 


[0052] 


Figui 


'e 36 is a hardware embodiment in which the inven- 




tion can operate. 



Detailed Description 



[0053] The present invention and the various features and ad- 
vantageous details thereof are explained more fully with 



reference to the non-limiting embodiments that are illus- 
trated in the accompanying drawings and detailed in the 
following description. It should be noted that the features 
illustrated in the drawings are not necessarily drawn to 
scale. Descriptions of well-known components and pro- 
cessing techniques are omitted so as to not unnecessarily 
obscure the present invention. The examples used herein 
are intended merely to facilitate an understanding of ways 
in which the invention may be practiced and to further en- 
able those of skill in the art to practice the invention. Ac- 
cordingly, the examples should not be construed as limit- 
ing the scope of the invention. 
[0054] First, the invention constructs a weighted Voronoi diagram 
(shown in Figures 1A and IB) for a particular electrical 
fault mechanism using known techniques. Denote this as 
Voronoi (fault). Electrical fault mechanisms are defined by 
the type of defect that may occur and one or more levels 
in the design involved in the fault. The Voronoi diagram 
that results maps points in the involved layout levels to 
their distance from edges of device shapes triggering the 
electrical fault (shown by the arrows in Figure IB). This 
Voronoi diagram partitions the layout into Voronoi cells 
110, or sets of contiguous points which share common 



fault-triggering edges 118 of a device shape 116. The 
boundaries 112 between cells 110 are known as Voronoi 
bisectors 112. The points where two or more Voronoi bi- 
sectors 112 meet are known as Voronoi vertices 114. 

[0055] a visualization of the mapping provided by Voronoi (fault) is 
that of a three-dimensional surface (e.g., see Figure 2). 
Imagine the layout design to be on the x-y plane and let 
the z-axis represent the mapped distance. The use of 
squares to model the shape of defects [1] results in 
Voronoi cells 110 that are piece-wise defined planar sur- 
faces, and bisectors 112 that are linear. This three- 
dimensional surface can be described by some planar, 
piece-wise continuous function over the positions and 
orientation of the bisectors 112 in the layout. 

[0056] Critical area represents the likelihood of a random defect 
and is a function over the surface the Voronoi diagram 
represents. The critical area within Voronoi (fault) is the sum 
of the critical areas contributions of its cells 110. Further- 
more, the critical area contribution of each cell can be 
computed as the sum of the critical area contribution of 
the bisectors bordering the cell. This is achieved using a 
trapezoidal decomposition technique. This is achieved us- 
ing a trapezoidal decomposition technique (as described 



in reference [l])The goal of yield improvement here is to 
reduce the critical area within the layout. While Voronoi 
(fault) allows us to compute the current critical area, it 
does not describe how critical area changes as a result of 
layout modification. In the following sections, a set of cost 
functions are presented which describe critical area in 
terms of variables representing the positions of edges 118 
in the design. This enables one to compute critical area as 
the edges 118 in the design change. This in turn enables 
one to optimize the geometry of the edges 118 in the de- 
sign to reduce critical area and, thus, improve yield. 
[0057] |_ e t k represent an edge 118 in the design as shown in 
Figure 2. 

Let x K = (x K ,y K ,0) represent a point (x K ,y K ) on K. 
Let g K ^(gfo.'gKy) represent the gradient of K. 

[0058] The gradient of K is the direction of increasing distance 

from K. The gradient is perpendicular to the orientation of 
K in the L-infinity distance metric, for all possible orienta- 
tions. Its magnitude is arbitrarily but consistently chosen. 

[0059] For the rectilinear boundary around the entire layout de- 
sign, we define special gradients as follows. 



^horizontal boundary ^(0,1) vertical boundary — ^(1,0) 



Let n K = (n Kx , n Ky ,n Kz ) represent the normal vector of the plane in three-dimensions 
representing distance away from K. For the L-infinity distance metric, 




Sk & {^(1,0)^(0,1) 
Sk = ^(1,0) 
8k = ^(o,i) 



(1,0,0) 
(0,1,0) 



[0060] | n other words, this portion of the disclosure defines the 
normal vector that is based on the orientation (encoded in 
the gradient) of the edges in the design. It shall be used 
to represent the mapping described by Voronoi cells in a 
way that is independent of the geometry comprising the 
Voronoi diagram itself. This can then be used to formulate 
cost functions to calculate the critical area contribution of 
the various edges of shapes within the integrated circuit 
design subject to modification and avoids the need to re- 
construct the Voronoi diagram under these modifications. 

[0061] As shown in Figure 3, Vertex (A, B, C) 114 represents the 
coordinate of a Voronoi vertex defined by edges A, B, and 
C of device shapes 116 in the design. This coordinate can 
be defined based on, but independent of, Voronoi vertices 
in the Voronoi diagram. Note that the devices, Voronoi bi- 



sectors, Voronoi vertices, etc. are not explicitly identified 
in all the drawings, so as to simplify the drawings and di- 
rect the reader's attention to the salient portions of the 
invention. 

Vertex (A, B, C) = ^ *^X« B * n c ) + (x B • n B ){n c x n A ) + {x c • n c ){n A x n B ) 

Det[n A n B n c ] 

[0062] • represents the dot-product operator 
[0063] x represents the cross-product operator 

[0064] | n other words, the coordinate of the vertex is a function 
of variables representing the positions and orientations 
(encoded as normals) associated with the involved edges 
in the design. This formulation is that of the intersection 
of the planes (as mentioned earlier) defined by the in- 
volved edges in the design. 

[0065] As shown in Figure 4, Bisector (A, B, C, D) 112 represents 
a Voronoi bisector whose geometry is defined by edges A, 
B, C, D of device shapes 116 in the design. These edges 
are stated in clockwise order around the bisector, such 
that B and D are on either side of the bisector. By defini- 
tion, Bisector (A, B, C, D) it is a line segment between ver- 
tices 114 Vertex (A, B, D) and Vertex (B, C, DJ.Let Contribution 
(Bisector (A, B, C, D) represent the critical area contribu- 



tion associated with a Voronoi bisector defined by edges 
A, B, C, and D in the design. Using squares to model the 
shape of defects, applying a trapezoidal decomposition 
technique (as described in Papadopoulou, E. and Lee, D.T., 
"Critical area computation via Voronoi diagrams, "Com- 
puter-Aided Design of Integrated Circuits and Systems, IEEE Trans- 
actions on , Vol. 18, No. 4, pp 463-474, April 1999.), and 
applying the equations derived above, critical area contri- 
bution is given by 

Contribution^Bisector (A, B, C, Z))) = 
k( NormalXY (B) - NormalXY(D)) x 



Vertex(C,B,D) xy -Vertex(A,B,D) xy Vertex {C,B,D) i 
Vertex(C,B,D) z -Vertex(A,B,D) z * Vertex (A, B,D) z / 



where k is a constant and 



NormalXY = 



'K*y 
"Kxy 



and for any vector f = (x,y, z) , 

%=( x >y) 



r=z 



[0066] | n other words, the critical area contribution of a Voronoi 
bisector is a function of the orientation (encoded as nor- 
mals) and positions of edges in the design. The critical 
area contribution of each Voronoi bisector is proportional 
to the cross product of two vectors the x-y difference 
vector between the normal vectors of the planes which 
meet at the bisector, and the x-y difference vector repre- 



senting the length of the bisector. The factor of propor- 
tionality is some constant times the logarithm of the z- 
coordinates of the vertices of the bisector divided by the 
difference in the z-coordinates of the vertices of the bi- 
sector. 

[0067] This equation is valid until the bisector or an adjacent bi- 
sector collapses to zero-length, resulting in a topological 
change in the Voronoi diagram. Having expressed the 
critical area contribution of a Voronoi bisector as a func- 
tion of edges in the design enables the invention to pre- 
dict change in critical areas as a result of layout modifica- 
tion. Based on the properties of the Voronoi diagram, the 
following observations can be made. Under continuous 
layout modification, the vertices in the Voronoi diagram 
shift positions predictably; the three-dimensional surface 
the Voronoi diagram represents changes in a continuous 
manner; the critical area of the layout changes continu- 
ously since it is a continuous function over this three di- 
mensional surface; and when such motion causes a 
Voronoi bisector to collapse to zero-length, its critical 
area contribution converges to zero: 

lim Contributioni Bisector (A, B, C, D)) = 0 

Vertex(C,B,D)-+Vertex{A,R,D) V V ') 



[0068] a further observation is that the point at which a bisector 
collapses is a function of the involved edges in the design. 
In fact, it is the point at which the coordinates of the two 
vertices defining the bisector converge. Therefore, one 
can predict the state in which a particular critical area 
contribution equation is valid under continuous modifica- 
tion of edges in the design. Furthermore, after collapse, 
an expansion may occur when modification continues, re- 
sulting in the emergence of a new bisector with a critical 
area contribution. This transformation is predictable and 
stated as follows (see Figure 5). 

Bisector (A,B,C,D) FWc^H^wfro ) {D ^ c) 

Bisector (A,B,C,D) I^-W^b^o > ^ g ^ ^ 

Bisector (A,B,C,D) I^^HW'^H ) B ^ ^ Q ^ 

Bisector (A, B, C, D) t^B^^M > Bisector {A, B, G, D) 

Bisector(A,B,C,D) t^^^M^ , Bisector (A, B^D) 

[0069] As illustrated, when Bisector (A, B, C, D) or a bisector ad- 
jacent to Bisector (A, B, C, D) in the Voronoi diagram col- 
lapses, a new bisector replaces Bisector (A, B, C, D). For 
example, suppose Vertex (B, C, D) and Vertex (H, C, D) are 
the first to converge. After convergence, Bisector (A, B, C, 
D) is replaced with Bisector (A, B, H, D) and therefore, 
Contribution (Bisector (A, B, C, D)) is replaced by Contri- 



bution (Bisector (A, B, H, D). 

[0070] N 0 te that the identification of the device shapes, etc. has 
been intentionally omitted from Figure 5 to more clearly 
illustrate the invention, as mentioned above. The replace- 
ment of one vertex for another is shown in the following 
transformations. 

[0071] For a Bisector (A, B, C, D) shown in Figure 5 the critical 
area cost function is stated as follows. 



Cost (Bisector (A,B,C, D)) = 



Contribution(Bisector (A,B,C,D)) Space I 
eContribution(Bisector (A, B, C, D)) Space II 



[0072] Contribution accurately describes the critical area contri- 
bution of the Voronoi bisector within Space I. This is the 
space in which the edges in the design start. Space II is 
the space of variables in which Contribution no longer ac- 
curately describes the critical area contribution of the 
Voronoi bisector due to a collapse of the bisector or of an 
adjacent bisector. These two spaces are separated by a 
multi-dimensional surface defined by 



Vertex (C, B, D) - Vertex(A, B, D)\\ = 0 
\\Vertex(E,B,A)-Vertex(D,B,A)\\ = 0 
\Vertex(F, A, D) - Vertex(B, A, D)\ = 0 
\\Vertex(D,B,C)-Vertex(G,B,C)\\ = 0 
\Vertex(B,C,D)-Vertex(H,C,D)\\ = 0 



[0073] where E, F, G, and H are those edges in the design relating 
to bisectors adjacent to Bisector (A, B, C, D) as in the ear- 
lier discussion. In the case where one accounts for topo- 
logical change through the transformations previously de- 
scribed, the value of e should be set to zero. Otherwise, 
the invention chooses a value e with which to penalize the 
cost function within Space II. By this technique, it confines 
the cost function to variables associated with only four 
edges at a time within the design. 

[0074] | n other words, instead of statically calculating critical 

area contribution from a Voronoi diagram, the invention 
utilizes a plane normal to represent each Voronoi cell de- 
fined by various device edges in the integrated circuit de- 
sign. Using these normals to calculate the Voronoi bisec- 
tors allows the invention to develop a cost function in 
terms of variables associated with the positions of device 
edges, independent of the Voronoi diagram itself. There- 



fore, when the positions of the device edges change, the 
values of each of the vertex coordinates will also change 
in accordance, thereby altering the critical area contribu- 
tion. By presenting the critical area contribution as a cost 
function, the invention allows an optimization process to 
reduce the critical area contribution cost. 

[0075] The invention associates this cost function with every 
Voronoi bisector Bisector (A, B, C, D) in the diagram. To 
extend the range of predictability, one can also associate 
cost functions with the possible transformations of Bisec- 
tor (A, B, C, D) described by the previous section. 

[0076] The critical area minimization objective is, therefore, 

Objective : Minimize 

X fault X Bisector B e Voronol(fau!t) X Bisector(A,B,C,D) flu transformations of B ^° St {^ SeCtor {A C, £>)) 

[0077] | n other words, the invention minimizes the sum of all 
costs from all bisectors, and transformations thereof 
(subject to constraints such as topological and ground- 
rule constraints), for every fault mechanism in the inte- 
grated circuit design. 

[0078] a non-linear, multi-dimensional optimization algorithm 
may be used to solve the foregoing objective for the posi- 
tions of all movable edges in the layout. Such algorithms 
are NP-complete or stochastic. Fortunately, there are a 



few simplifications that can be made to the cost function 
that would enable the use of a linear optimization algo- 
rithm, which is considerably more efficient. 

[0079] |_ e t a, b, c, d, be some constants. Let w represent the sum 
or difference of two variables; the set of variables being 
the x and y coordinate of all movable edges in the layout. 
Without loss of generality, apply the following simplifica- 
tions. The example used herein assume orthogonal data, 
however, one ordinarily skilled would understand that any 
arbitrary angles could be used. 

[0080] |_ Restrict the set of variables to represent a single dimen- 
sional movement (other dimensions may be handled using 
various techniques, for example, constraining them to 
move with an adjacent orthogonal edge.) Then 



[0081] The critical area contribution of a Voronoi bisector simpli- 
fies to 



^€{(-l,0),(l,0),(0,-l),(0,l),/? (10> ^ {01) } 

n K e {(1,0,1),(-1,0,1),(0,1,1>,(0,-1,1>,(1,0,0),{0,1,0)} 



Contribution [Bisector (A,B,C,D))e\a^+b, a\n^- + b 




[° 082 ] II. Ignore the effect of adjacent Voronoi bisectors collaps- 



ing and consider only the Voronoi bisector itself collaps- 
ing. Then the condition at which Contribution is no longer 
accurate can be expressed by 

cw, + dw 2 = 0 

[° 083 ] III. Solve the objective in two passes. In the first pass, 
consider only those variables which represent x- 
coordinates. In the second pass, consider only those vari- 
ables which represent y-coordinates. Then 

Contribution(Bisector (A, B, C, D)) e jmv + b, — + Z>, a In w + b, a In vt>, + b In w 2 j 

[0084] iv. Approximate Contribution using a piecewise-linear 
function. 

[0085] Actions I and II result in a cost function over sums or dif- 
ferences of variables. Actions III and IV result in cost func- 
tions that are convex and piecewise-linear. An efficient 
optimization algorithm such as graph-based simplex can 
then be applied. 

[0086] | n other words, allowing edges to move along a single axis 
(x-axis or y-axis) simplifies the cost function and enables 
it to be used within a linear optimization algorithm which 
is considerably more efficient that a general optimization 
algorithm. 



[0087] There are a number of topologies that may exist within 
the Voronoi diagram which are not handled by the above 
cost function. These fall into one of two categories ambi- 
guities and degeneracies. Ambiguities shall now be ad- 
dressed. 

[0088] An ambiguity is present when the plane normals on either 
side of a Voronoi bisector are parallel to each other as 
shown in Figure 6A. For example, 

Bisector ({P,Q} ,B,C,D) may have the property that g B = g D and, consequently, 
n B =ii D . Functions Vertex ({P,Q},B,D) and Vertex (C,B,D) are only applicable when 
P, Q, B, D and C, B, D have non-parallel plane normals. For orthogonal data, 
modification to bisectors B and D may result in one of two scenarios shown in Figures 
6B and 6C where a Voronoi cell emerges at the bisector. 

[0089] jo account for such changes, the following decomposition 
can be used 



Bisector ({P, Q) , B, C, D) -> < 
where g p *g Q ,g p ± g B ,andg Q 1 g B 



Bisector (Q, P, C, D) , Bisector (D, P, B, C) , 
Bisector (P, B, C, Q) , Bisector (D, Q, B, C) 



[0090] | n a scenario such as this, where one accounts for topo- 
logical change, the value e in the cost functions associated 
with bisectors resulting from the decomposition should be 
set to zero. 

[0091] | n other words, when a bisector is encountered with ambi- 



guity, the invention decomposes the bisector to its poten- 
tial forms during layout modification, resulting in unam- 
biguous bisectors. This rule may be applied iteratively un- 
til no ambiguities are present. A cost function is then de- 
rived for every Voronoi bisector at the end of the decom- 
position. Under layout modification the cells 110 that 
emerge are initially of zero area and contribute no critical 
area to that of the overall layout. If one scenario becomes 
dominant, its associated cost functions will assert them- 
selves within the space of predictability while those of the 
other scenarios are suppressed as a consequence of being 
outside the space of predictability. The overall effect is an 
accurate modeling of how critical area will change under 
layout modification in cases of ambiguity. 

[0092] As mentioned above, there are a number of topologies 
that may exist within the Voronoi diagram which are not 
handled by the cost function and these fall into one of two 
categories ambiguities and degeneracies. Degeneracies 
shall now be addressed. 

[0093] a degeneracy is present when more than three bisectors 
converge at a single Voronoi vertex, as shown in Figure 
7A. Modification to the layout may result in an emergence 
of one or more bisectors at this vertex, as shown in Fig- 



ures 7B-7G. In other words, the invention decomposes 
these vertices into all possible scenarios that may result 
under layout modification, resulting in additional vertices 
and bisectors. Figures 7B-7G illustrate the different possi- 
bilities that can occur from Figure 7A under layout modifi- 
cation. Thus, for every subset of four edges in the design 
defining the Voronoi vertex, the invention applies the fol- 
lowing decomposition. 

Vertex({A,B,C,D}) -+ {Bisector (A, B,C,D), Bisector (B,C,D, A)} 

[0094] Additionally, the invention derives a cost function for ev- 
ery Voronoi bisector resultant from this decomposition. 
The Voronoi bisectors in this decomposition initially have 
zero-length and their critical area contribution is zero. As 
the bisectors grow in length, their cost functions move 
into the space of predictability and dominate those asso- 
ciated with Voronoi bisectors that do no emerge. In a sce- 
nario such as this, where one accounts for topological 
change, the value e in the cost functions associated with 
bisectors resulting from the decomposition should be set 
to zero. The overall effect is an accurate modeling of how 
critical area will change under layout modification in cases 
of degeneracy. 



[0095] Figure 8 illustrates the flow for yield improvement using 
the critical area minimization objectives. More specifically, 
in item 800, the invention inputs an initial layout. For ev- 
ery applicable fault mechanism (parameterized by type of 
defect and the involved layout levels), the invention builds 
layout variables for movement in one axial direction (as 
shown in item 802). More specifically, the invention asso- 
ciates variables to the positional components of the edges 
that lie in a first direction in the initial integrated circuit 
design. Then, in item 804, the invention constructs 
Voronoi diagrams. The creation of these Voronoi diagrams 
occurs incrementally as processing moves across the inte- 
grated circuit design. As each Voronoi bisector is final- 
ized, processing moves to item 806 which decomposes 
the bisector (eliminates ambiguity or degeneracy as dis- 
cussed above). 

[0096] Once each Voronoi bisector is decomposed, the invention 
associates a cost function in terms of critical area contri- 
butions of the Voronoi bisectors in item 808. As shown 
above, this cost function is a function of the positions and 
orientations (encoded as normals) of the edges, and the 
critical area contributions comprise a measure of electrical 
fault characteristics of the areas between the edges. 



[0097] | n items 810, the invention applies a linear optimization 
algorithm to all the cost functions that are collected from 
item 808 and processing optionally returns to item 800 to 
optimize for another axis (orthogonal, or any other useful 
direction (45 degrees, 60 degrees, etc.)). Thus, the inven- 
tion optimizes the positions and length of the edges to 
reduce critical area contribution cost in the first direction 
across the integrated circuit design to produce a revised 
integrated circuit design and then, the invention repeats 
this process with the revised integrated circuit design in a 
different direction. 

[0098] | n a standard run, optimization is performed twice, once 
along each of the x and y axes. For further yield improve- 
ment, more iterations of this flow can be executed. This 
flow may be combined with the generation of other cost 
functions and objectives such as topological or ground- 
rule constraints and objectives. When combining multiple 
objectives, each should be weighted appropriately relative 
to each other. 

[0099] Figures 9A-32C illustrate some possible Voronoi bisector 
configurations. More specifically, the "A" Figures (of Fig- 
ures 9A-32C) illustrate the specific Voronoi bisector con- 
figuration and the appropriate cost functions for such 



Voronoi bisector configuration. In the "A" Figures, the ar- 
rows represent the direction of gradient. The "B" Figures 
(of Figures 9A-32C) illustrate the graphical results of the 
cost function in one direction and the "C" Figures (of Fig- 
ures 9A-32C) illustrate the graphical results of the cost 
function in a direction orthogonal to the first direction. 
The illustrated cost functions also apply to configurations 
that are rotationally or reflectively symmetric to the drawn 
configuration. 

[0100] Figures 33-35 represent Voronoi bisector configurations 
which do not contribute a cost function to the optimiza- 
tion algorithm because modifications of the involved vari- 
ables do not change their critical area contributions. For 
this reason, they are ignored and not included within the 
total of all cost functions (item 808). 

[0101] a representative hardware environment for practicing the 
present invention is depicted in Figure 36, which illus- 
trates atypical hardware configuration of an information 
handling/computer system in accordance with the subject 
invention, having at least one processor or central pro- 
cessing unit (CPU) 10. CPUs 10 are interconnected via sys- 
tem bus 12 to random access memory (RAM) 14, read- 
only memory (ROM) 16, an input/output (I/O) adapter 18 



for connecting peripheral devices, such as disk units 11 
and tape drives 13, to bus 12, user interface adapter 19 
for connecting keyboard 15, mouse 17, speaker 103, mi- 
crophone 104, and/or other user interface devices such as 
touch screen device (not shown) to bus 12, communica- 
tion adapter 105 for connecting the information handling 
system to a data processing network, and display adapter 
101 for connecting bus 12 to display device 102. A pro- 
gram storage device readable by the disk or tape units, is 
used to load the instructions which operate the invention 
also loaded onto the computer system. 

[0102] Therefore, as shown above, the invention presents a 

method, system, program storage device, etc. that calcu- 
lates and reduces critical area in an integrated circuit de- 
sign. More specifically, the invention begins with an inte- 
grated circuit design as its initial input. It associates vari- 
ables with the positions of edges in the integrated circuit 
design. It encodes the orientation of edges in the inte- 
grated circuit design as normals. 

[0103] The invention defines cost functions of these variables in 
terms of the critical area contributions within the spacings 
for the edges associated with these variables. The critical 
area contributions comprise a measure of the electrical 



characteristics of the integrated circuit for a particular 
type of fault. 

[0104] The invention associates cost functions of these variables 
with particular spacings among particular edges in the 
design. This process uses a three-dimensional represen- 
tation of the Voronoi diagram constructed for the inte- 
grated circuit design for a particular type of electrical 
fault. In this process, the invention defines Voronoi cells 
as a mapping of points in the spacing between edges in 
the design to the minimum size of defects that trigger an 
electrical fault between them. The set of points where 
Voronoi cells meet define Voronoi bisectors. The end- 
points of these bisectors form Voronoi vertices. 

[0105] The invention represents Voronoi cells, bisectors, and ver- 
tices using mathematical expressions, independent of the 
Voronoi diagram itself, over the variables associated with 
the edges of the integrated circuit design involved in the 
spacing. Based on these entities, the invention formulates 
cost functions which describe the critical area contribu- 
tions of the spacing under continuous modification of the 
design. 

[0106] The invention uses these cost functions in an optimization 
algorithm to modify the positions and length of the edges 



in the integrated circuit design to reduce critical area con- 
tribution cost. For improved efficiency with a linear opti- 
mization algorithm, the invention applies the optimization 
for variables in one direction across the integrated circuit 
design to produce a revised integrated circuit design. The 
invention optionally repeats this process using the revised 
integrated circuit design in another direction to further 
reduce critical area cost. 

[0107] As we enable ourselves to build smaller and denser inte- 
grated circuits, their sensitivity to the occurrence of spot 
defects inherent in the manufacturing process becomes of 
greater importance. Such defects cause electrical faults 
within the circuit, contribute to yield loss, and ultimately 
result in lost resources. This invention addresses the 
problem by providing an automated means of optimizing 
a circuit design to reduce this sensitivity. Furthermore, the 
process described herein can be combined with other op- 
timization objectives to achieve a balance of costs in the 
modified design. 

[0108] while the invention has been described in terms of pre- 
ferred embodiments, those skilled in the art will recognize 
that the invention can be practiced with modification 
within the spirit and scope of the appended claims. 



